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Method of determining filter coefficients from line spectral frequencies 




\ 

The present invention relates to a method of determining filter coefficients 
firom Line Spectral Frequencies comprising recomputing P(z) and Q(z) polynomials and 

\ 

comprising calculating the co^ coefficients. 

The coding of speech signals is used particularly in the field of mobile 
communications since the coded speech signal can be transmitted in a manner in which the 
redundancy commonly experienced in human speech is reduced. Linear Predictive Coding 
® (LPC) is a known technique normally used in speech coding and in which the correlation of 

%j the speech signal is removed by means of a filter. The filter is best described by way of one 

^ of a different set of parameters, and one important set of which comprises LSFs. 

yi 10 An accurate representation of the filter is an important requirement since such 

^ information is transmitted with the speech signal for subsequent reconstruction of the speech 

_2 signal at a signal-receiving unit. 

fy The advantages of representing LPC filter coefficients in the form of LSFs 

2 have been well-documented since the inception of this concept in 1975. As is well known, 

15 the representation of an inverse LPC filter A(z) in the form of LSFs is derived from the 
representation of A{z) by its set of zeros in the z-plane. Insofar as the function A(z) 
represents an all-zero filter, it can be fully and accurately described by way of reference to its 
corresponding set of zeros. 

Computation of the LSFs commences with the decomposition of the 
20 polynomial (z) of order m into two inverse polynomial functions P{z) and Q{z) . For 
confirmation, the polynomial A^ (z) and the two inverse polynomials appear as follows: 
A^ (z) = 1 + a,z-' + a^z-^ + ... + a^z"" 
and 

i>(z) = ^„(z) + z-('"*'>^„(z-') 
25 Qiz) = A„(z)-z-^-*'U„iz-') 

The polynomials P(z) and Q{z) each have (m + 1) zeros and exhibit various 
important characteristics. In particular: 
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all zeros of P{z) and g(z) are found on the unit circle in the z-plane; 
the zeros of P(z) and Q(z) are interlaced on the unit circle and the zeros do 
not overlap; and 

the minimum phase property of (z) is easily preserved when the zeros of 
P(z) and Q{z) are quantised. 

Analysis of the above confirms that z = -1 and z = +1 is always zero with the 
functions P(z) and Q(z) and since these zeros do not contain any information relating to the 

LPC filter, they can simply be removed from P(z) and Q(z) by dividing by (1 + z~* ) and 

(1-z-^). 

Such revised functions can be represented when m is even as follows: 
(1 + z-') (1-^'') 

and when m is odd as: 



P' (z) = P(z) and 2' (z) = - 

^ ^ ^ ^ (l-z-')(l + z-*) 

The advantageous properties of functions P(z) and Q(z) as noted above are 
also valid for P\z) and g' (z) . Since the coefficients of P'(z) and Q\z) comprise real 
numbers, the zeros form complex conjugate pairs such that the search for zeros only has to be 
conducted on the upper half of the unit circle, i.e. where 0<o)<7r , 

It generally proves inconvenient to compute complex zeros, particularly by 
way of computerised numerical analysis methods, and so the functions P' (z) and Q' (z) are 
transformed to functions P"(z) and Q"(z) with real zeros. Also, the functions F(z) and 
Q' (z) always have an even order and, since they are symmetrical, the functions can be re- 
written with real zeros to the following manner: 

(CO) = 2^ " cos((m^ - 0^) 
g'(6>) = 2j;^,"cos((/w^ -i)co) 

(=0 

where Po''=U Pxx..m,-x'^ PxX..m,-X ^ P., "=|Pm/.^o"= U ^l,2...m,-l"=^1.2...m,-l'. 

q '':= — q \ and where m „ is equal to the number of zeros of P* (z) on the upper half of 

q 2 *f 
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the unit circle and where is equal to the number of zeros of Q' (z) on the upper half of the 
unit circle. 

When seeking the zeros of these functions, advantage can be taken from the 
form of the representations for P"(z) and due to the fact that the number of zeros to 

be located is already known. 

Importantly, and of particular relevance to the present invention, once the 
LSFs have been identified and employed as required, the recomputation of the LPC filter 
coefficients from the LSFs is required. While this stage represents a much less 
computationally intensive calculation than the computation of the LSFs from the filter 
coefficients as discussed above, problems and disadvantageous limitations are experienced. 
In particular, the values of intermediate coefficients can be disadvantageously high and this 
can lead to numerical problems even when employing floating point representations. 

Reconibuting LPC filter coefficients a, from LSFs is much less 
computationally intensive than computing the LSFs from the filter coefficients. Each LSF 
0). , / = 0,1,..., m - 1 contributes to a quadratic factor of the form, 1-2 cos(<s)^ )z"^ + z"^ . The 
polynomials P'(z) and Q\z) are formed by multiplying these factors using the LSFs that 
come from the corresponcKng polynomial: 

mp-\ ] 

/>'(z) = n(\-2cos(o2,)z-' +2-') 

1=0 



Q' (z) = fl (1 \ 2 cosia,,,,, )z-'+z-') 

1 = 1 



The polynomials P(z) and g(z) are computed by multiplying P'(z) and Q\z) with the 
extra zeros at z = -1 and z = +1\. Finally, the filter coefficients are computed by using the 
following equation: 

2 

which defines the relationship between the polynomial (z) and the two inverse 

polynomials discussed earlier. 

Thus, when recomputing the P(z) and g(z) polynomials one can use the 

above equations for P'(z) and Q'iz) and add the extra zeros. Thus, for m is even: 

mp-\ 

P{z) = (1 + z ) fl (1 - 2 cos(<W2, )z-' + ) 

1=0 
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Q(z) = (1 - )n (1 - 2 cosC^y^,., + ^'^) 
and for m is odd: 

P(z) = ]^(l-2cos(fi>2,)z-' +2'') 

1=0 

e(z) = (1 - )(1 + )n (1 - 2 cosK,,, )z-' + ) 

Since ft?, coefficients are ordered in increasing frequency the contributions of the first 
cos(cOf ) coefficients are positive and the last coefficients are negative. This introduces an 
undesirable growth of the intermediate coefficient values while conducting the polynomial 
multiplications (1 - 2cos(a>2/+i )^ '^ + z"^) . With an increasing order of aw, such problems are 
amplified. To illustrate this, an example polynomial Q(z) = 1 - z"^^ has been taken; it being 
noted that that m = 2N , Such a polynomial provides for a system with 2N equidistant zeros 
on the unit circle. Since this is only a very simple example, it should of course be appreciated 
that in reality the growth of the intermediate coefficients can be much larger. Therefore a 
different strategy must be used. It has been found that for high order polynomials m=60 (or 
N=30) a double precision floating point representation will not be sufficient. A different 
method might involve searching for the best possible combination of co^ and which exhibits 
the least amount of intermediate coefficient growth. However, due to the large number of 
possible combinations, this is most unlikely to be feasible and it also means that the optimal 
combination will never be found. 

The invention seeks to provide for a method of determining filter coefficients 
having advantages over known such methods. 

In accordance with the present invention there is provided a method of 
determining filter coefficients from Line Spectral Frequencies as noted above and 
characterised by the steps of addressing the polynomials in a series and reducing the number 
of polynomials in co^ in the said series by combining the polynomials in co^ two by two in a 
manner so as to arrive at two polynomials in co^ and determining the product of the said two 
polynomials. 

The invention serves to combine cd^ in such a way that hardly any signal 
growth occurs and proves particularly advantageous since the use of an increasing index i 
would not seem to offer a good solution. 



m 
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Using the method of the present invention with the example polynomial noted 
above, i.e. Q(z) = 1- z"^^ , the intemiediate coefficients are never larger than 2. In practice 
only a limited amount of intermediate coefficient grov^h occurs. Advantageously, the 
invention need not comprise a particularly complex method. In general it requires only a 
different indexing and can advantageously deliver almost optimal results. For P(z) the same 
procedure could be used only if m is even then P(z) has a root at z = -1 . If w is odd P(z) 
does not have any additional roots, so no additional roots need to be added. 

Finally the relationship noted above is employed to compute the coefficients 
of A(z) from P(z) and Qiz) . 




^^^^T^ The invention is described further hereinafter, by way of example only, with 
reference to the accompanying drawing which is a graphical representation of the 
intermediate coefficient growth experienced in the prior art and in an example polynomial 



15 Q(z) = l-z 



-2N 




20 




it is assumed that the original polynomial is reconstructed by combining the 
zeros with increasing co^ . The maximum value of the largest coefficient during the 
recomputation pfbcedure is plotted in the accompanying drawing. Note that the Y axis is 
logarithmic. For l^ge order the intermediate values of some of the coefficients become 
very high. 

However, such problems are not encountered in a method according to the 
present invention. 

As an example, and for Q{z) vAth m is even, the following ordering of the 
polynomials is uset 

Vo[OK=l-^" 



Vo[l] =\l - 2cosft?,z * -f z 
v^lll^iSr 2 cos CO ^z'^ +z 
30 Vo[mJ = l\2cos6;2.„^_, 

If m is odd then the terms fo? Q{z) are 



-2 



Z + Z 
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Vo[0] = l-z-> 

Vo[l] = l-2cOS6>iZ"' 4-Z~^ 
Vo['wj = l-2COS6;2*^^_,Z-' +Z-' 

The next step is to combine the polynomials Vq [/] . The strategy will be 
demonstrated with an example with m=J2 and mg=6. The original seven polynomials are: 
vofOJ. vo[l], vo[2], vo[3], vo[4], vo[5], and vo[6]. 

In the first step the polynomials are combined two by two. Polynomial /' is 
combined with polynomial [mq-i], this gives four intermediate polynomials v,///* 
O 10 v,[0] = vJO].Vo[6] 

v,[l] = Vo[l]-Vo[5] 
v.[2] = Vo[2].Vo[4] 
v.[3] = vJ3] 

These four polynomials are combined in the same way, leading to two new 

\ 15 polynomials v^/y/ * 

v,[0] = v,[0].vi[3] 

v,[l] = vJl].vJ2] 
^^NTY^^ — yr The proclVict V2[0],V2[1] gives the final result: 
^ / V3[0] = V2)^].v,[l] 
20 The procedure can formally be described by the following pseudo program: 
if m is even 
begin 
mq=m/2 
mc=mq+l 
25 end 
else 
begin 

mq=(m-l)/2 
mc=mq-^2 
30 end 

np=mc 
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i=mc»I 



/*arithmetic shift right*/ 



k=0 

while (i>0) 
begin 

n^O 

while (n<i) 
begin 

Vk^i[n]=Vk[n], Vk[np-n-l] 
end 

if Hp is odd then 
begin 

Vk+i[n]=Vk[n] 

n=n~^l 

end 

np=n 

k=k+l 

i=n> >L /*arithmetic shift right */ 



Using this method with the example polynomial Q{z) = 1 - z" the 



intermediate coefficients are never larger than 2. In practice only a limited amoimt of 
intermediate coefficient growth occurs. This is not a very complex method (actually it only 
uses a different indexing) and delivers almost optimal results. For P{z) the same procedure 
will be used only if m is even then P{z) has a root at z = -1 . If w is odd P{z) does not have 
any additional roots, so no additional roots need to be added. The last step consists of using 
the equation 



to compute the coefficients of A{z) from P(z) and Q{z) . 



end 



